Preprint typeset in JHEP style. - HYPER VERSION 



HUTP-02/A029 
hep-th/0207096 



Chern-Simons theory, matrix integrals, 

o and perturbative three-manifold invariants 

o 

> 

O 



Marcos Marino 

Jefferson Physical Laboratory, 
Harvard University, 
Cambridge, MA 02138, USA 



> 

\^ • marcos@born.harvard.edu 

o\ 
o 
l> 
o 

CN Abstract: The universal perturbative invariants of rational homology spheres can be extracted 

^ from the Chern-Simons partition function by combining perturbative and nonperturbative results. 

^ • We spell out the general procedure to compute these invariants, and we work out in detail the case 

• of Seifert spaces. By extending some previous results of Lawrence and Rozansky, the Chern-Simons 

D ' partition function with arbitrary simply-laced group for these spaces is written in terms of ma- 

, trix integrals. The analysis of the perturbative expansion amounts to the evaluation of averages 

' in a Gaussian ensemble of random matrices. As a result, explicit expressions for the universal 

^ perturbative invariants of Seifert homology spheres up to order five are presented. 
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1. Introduction 

Chern-Simons theory has been at the heart of the developments in three-manifold topology 
and knot theory for the last ten years. The partition function of Chern-Simons theory defines 
a topological invariant of three-manifolds, sometimes known as the Witten-Reshetikhin-Turaev 
invariant, that can be studied from many different points of view. In general, the invariant thus 
obtained contains information about the three-manifold itself but also about the gauge theory group 
that one uses to define the theory. 

However, from a perturbative point of view it is clear that one can extract numerical invariants 
of the three-manifold which are intrinsic to it and do not depend on the gauge group. This goes 
as follows: if we compute the partition function in perturbation theory, the contribution at a given 
order consists of a sum of terms associated to Feynman diagrams. Each term is the product of a 
group dependent factor (the group weight of the diagram) , and a factor involving multiple integrals 
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of the propagators over the three-manifold. This last factor does not depend on the gauge group one 
started with, and in this sense it is universal. Therefore, one can extract from perturbation theory 
an infinite series of invariants, the so-called universal perturbative invariants of three-manifolds. 

The idea of looking at the perturbative expansion of Chern-Simons theory in order to extract 
numerical invariants that "forget" about the gauge group was first implemented in the context of 
knot invariants, leading to the theory of Vassiliev invariants and the Kontsevich integral (see |^, |28[]). 
The perturbative approach to the study of the partition function of Chern-Simons theory has a long 
story, starting in This has been pursued from many points of view. On the one hand, the 

structure of the perturbative series has been analyzed in detail (see for example 0, ^ and |16| for a 
nice review) , leading to the graph homology of trivalent graphs as a systematic tool to organize the 
expansion. On the other hand, the asymptotic expansion of the nonperturbative results has also 
been studied PBI, 130, ETl, EDI, |3D|, although so far all the analysis have focused on theories with 



gauge group SU{2). Finally, a mathematically rigorous theory of universal perturbative invariants 
of three-manifolds has been constructed starting from the Kontsevich integral: the so-called LMO 



invariant ||31]] and its Aarhus version [g. 

The main goal of the present paper is to elaborate on the topological field theory approach 
to universal perturbative invariants. The point of view presented here is very similar to the one 
advocated in [0, ^ to extract Vassiliev invariants from Chern-Simons perturbation theory: first, 
one analyzes the structure of the perturbative series of an observable in the theory. This means 
in practical terms choosing a basis of independent group factors and compute its value for various 
gauge groups. In a second step, one computes the corresponding invariant nonperturbatively for 
those gauge groups, performs an asymptotic expansion, and extracts the universal invariants by 
comparing to the perturbative result. This program was applied successfully in [|], |[ to compute 
Vassiliev invariants of many knots. It turns out that, in the case of the Chern-Simons partition 
function, the first step is relatively easy, but the calculation of the partition function for arbitrary 
gauge groups in a way that is suitable for an asymptotic expansion turns out to be trickier, except 
in very simple cases. 

In this paper, some well-known results concerning the structure of the perturbative series are put 
together, and we carry out a detailed analysis up to order five. The focus is on a rather general class 
of rational homology spheres, Seifert spaces. The partition function of Chern-Simons theory with 
gauge group SU{2) on these spaces and its asymptotic expansion have been studied in |2^, p7| . 



The extension to higher rank gauge groups has also been considered HO, E^, but in forms that 



are not useful for a systematic perturbative expansion. In [^, Lawrence and Rozansky found a 



beautiful expression for the SU{2) partition function on Seifert spaces in terms of a sum of integrals 
and residues. It turns out that their result can be generalized to any simply-laced group and written 
in terms of integrals over the Cartan subalgebra of the gauge group (these kind of integrals already 
appeared in a related context in |36]). Interestingly, they are closely related to models of random 
matrices, and one can use matrix model technology to study the Chern-Simons partition function 
on these spaces. The resulting expressions can be expanded in series in a fairly systematic way, 
and by comparing the result with the general structure of the perturbative expansion, the universal 
perturbative invariants can be extracted. It should be mentioned that the full LMO invariant of 
Seifert spaces has been computed by Bar-Natan and Lawrence by using techniques from the 
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theory of the Aarhus integraL However, their result is rather imphcit and involves a complicated 
graphical calculus. 

This paper is organized as follows: in section 2, we review the computation of the Chern-Simons 
partition function starting from a surgery presentation. In section 3, we analyze in some detail the 
structure of the Chern-Simons perturbation series. In section 4 we compute the exact partition 
function of Seifert spaces for simply-laced gauge groups, generalizing the results of Lawrence and 
Rozansky, and we make the connection to matrix models. In section 5, we analyze the asymptotic 
expansion of the exact result, explain how to evaluate the matrix integrals, and present the results 
for universal perturbative invariants up to order five. In section 6, we comment on the possible 
relevance of these results to other physical contexts, and some avenues for future research are 
suggested. The Appendix collects the explicit expressions for the group factors and the matrix 
integrals, together with a summary of the properties of symmetric functions that are used in the 
paper. 



2. The partition function of Chern-Simons theory 

In this section we review some well-known results about the computation of the Chern-Simons 
partition function in terms of surgery presentations. An excellent summary, that we follow quite 



closely, is given in [36 



We consider Chern-Simons theory on a three-manifold M and for a simply-laced gauge group 
G, with action 

S{A) = ^ ! ^I(A^dA + -A^A^A\ (2.1) 



4vr Jm ^ 3 

where A is a G-connection on M. We will be interested in framed three-manifolds, i.e. a three- 
manifold together with a trivialization of the bundle TM © TM. As explained in ||^ , for every 
three-manifold there is a canonical choice of framing, and the different choices are labeled by an 
integer s G Z in such a way that s = corresponds to the canonical framing. Unless otherwise 
stated, we will always work in the canonical framing, and we will explain below how to incorporate 
this in the calculations, following |2^, |50 |. 



As shown by Witten in |Q, the partition function of Chern-Simons theory 

Zk{M) = I PAe^^cs(A)_ (2.2) 



defines an invariant of framed manifolds. There is a very nice procedure to evaluate ( p.2| ) in a 
combinatorial way which goes as follows. By Lickorish theorem (see for example ||32|), any three- 
manifold M can be obtained by surgery on a hnk C in S^. Let us denote by /Cj, i = 1, ■ ■ ■ ,L, 
the components of C. The surgery operation means that around each of the knots /Cj we take a 
tubular neighborhood Tub(/Cj) that we remove from S^. This tubular neighborhood is a solid torus 
with a contractible cycle and a noncontractible cycle jSi. We then glue the solid torus back after 
performing an SL(2, Z) transformation given by the matrix 

jj{p.,q.) = fPiri\ 1^2.3) 

\ Qi Si J 
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This means that the cycles PiOi + qi(3i and r^ai + Si/3i on the boundary of the complement of ICi are 
identified with the cycles aj, Pi in Tub(/Cj). 

This geometric description leads to the following prescription to compute the invariants in 
Chern-Simons theory. By canonical quantization, one associates a Hilbert space to any two- 
dimensional compact manifold that arises as the boundary of a three-manifold, so that the path- 
integral over a manifold with boundary gives a state in the corresponding Hilbert space. As it was 
shown in [Q, the states of the Hilbert space of Chern-Simons theory associated to the torus are 
in one to one correspondence with the integrable representations of the WZW model with gauge 
group G at level k. We will use the following notations in the following: r denotes the rank of G, 
and d its dimension, y denotes the dual Coxeter number. The fundamental weights will be denoted 
by Aj, and the simple roots by with i = 1, ■ ■ ■ ,r. The weight and root lattices of G are denoted 
by Aw and Ar, respectively. Finally, we put I = k + y. 

A representation given by a highest weight A is integrable if the weight p + A is in the fun- 
damental chamber JF; [p denotes as usual the Weyl vector, given by the sum of the fundamental 
weights). The fundamental chamber is given by A„//Ar modded out by the action of the Weyl 
group. For example, in SU{N) a weight p = X]I=iPi-^« -^i 



< I, and Pi > 0, i = 1, 



r. 



(2.4) 



i=l 



In the following, the basis of integrable representations will be labeled by the weights in Ti. 

In the case of simply-laced gauge groups, the Sl(2, Z) transformation given by [/(P''') has the 
following matrix elements in the above basis E^, : 



a/3 



\i sign(g)]l^+l 

7^exp 



Vol Aw 
VolAr 



^ ^ e{w) exp|^(pa2 - 2a{ln + w{P)) + s{ln + w(/?))^|. 



nGAr/(?Ar w6W 



(2.5) 



In this equation, |A_|_| denotes the number of positive roots of G, and the second sum is over the 
Weyl group W of G. ^{U^''^^) is the Rademacher function: 



p r 
q s 



p + s 



12s{p,q), 



(2.6) 



where s{p, q) is the Dedekind sum 



1 -nn 

71=1 



imp 



(2.7) 



With these data we can already present Witten's result for the Chern-Simons partition function 
of M. As before, suppose that M is obtained by surgery on a link £ in S^. Then, the partition 
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function of M is given by: 

Z(M,0 =e^^- Z^,^...^Ml^iT^---l^tT^- (2.8) 

In this equation, ... q^(£) is the invariant of the hnk C with representation ai — p attached to 
its z-th component (recall that the weights in J-'i are of the form p + A). The phase factor e**^'' is a 
framing correction that guarantees that the resulting invariant is in the canonical framing for the 
three-manifold M. Its explicit expression is: 

^fr = ^(j2^iU^''^''^)-SaiC)), (2.9) 

where o"(£) is the signature of the linking matrix of C. 



3. Chern-Simons perturbation theory 



The expression ( p.8|) gives the nonperturbative result for the partition function of M, and allows an 
explicit evaluation for many three-manifolds for any gauge group G and level k. However, from the 
point of view of Chern-Simons perturbation theory, the partition function can be also understood as 
an asymptotic series in whose coefficients can be computed by evaluating Feynman diagrams. 
In this section we review some known facts about the perturbative expansion of Chern-Simons 
theory and we state our strategy to compute the universal perturbative invariants. 

We are interested in the perturbative evaluation of the partition function ( p.2|) . Let us assume 
(as we will do in this paper) that M is a rational homology sphere. The classical solutions of the 
Chern-Simons action are just fiat connections on M, and for a rational homology sphere these are 
a finite set of points. Therefore, in the perturbative evaluation one expresses Zk{M) as a sum of 
terms associated to stationary points: 

Z,{M) = J2z^\m), (3.1) 

c 

where c labels the different fiat connections A^'^'' on M. Each of the terms in this sum has a 
perturbative expansion as an asymptotic series in The structure of the perturbative series was 
analyzed in various papers ^ ^] and is given by the following expression: 



(M) = Z;i„p(M). exp J2 Sf'^' • (3.2) 



In this equation, x is the effective expansion parameter: 



x=—. (3.3) 
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The one- loop correction z[^}^^^^(M) was first analyzed in |44], and has been studied in big detail 
since then. It has the form, 



{2nx 



i(dimHO_dim//l) 



T. 



R I 



(3.4) 



where H^'^ are the cohomology groups with values in the Lie algebra of G associated to the flat 
connection t]^^ is the Reidemeister- Ray-Singer torsion of A^''\ is the isotropy group of 
A(^), and (y9 is a certain phase. More details about the structure of this term can be found in 

2o|, ^ g. 

Our main object of concern in this paper are the terms in the exponential of (|3.2|) corresponding 



to the trivial connection, which we will simply denote by S^. In order to make a precise statement 
about the structure of these terms, we have to explain in some detail what is the appropriate set of 
diagrams we want to consider. In principle, in order to compute S"^ we just have to consider all the 
connected bubble diagrams with d loops. To each of these diagrams we will associate a group factor 
times a Feynman integral. However, not all these diagrams are independent, since the underlying 
Lie algebra structure imposes the Jacobi identity: 



^ ^j i^fabefedc ~l~ fdaefebc ~l~ facefedb) 0. 



(3.5) 



This leads to the diagram relation known as IHX relation. Also, antisymmetry of fabc leads to the 
so-called AS relation (see for example ^ |TB|, ^). The existence of these relations between 
diagrams suggests to define an equivalence relation in the space of connected trivalent graphs by 
quotienting by the IHX and the AS relations, and this gives the so-called graph homology. The space 
of homology classes of connected diagrams will be denoted by ^(0)'^°'^°. This space is graded by 
half the number of vertices, and this number gives the degree of the graph. The space of homology 
classes of graphs at degree i is then denoted by A{^)f'^^. For every i, this is a finite-dimensional 
vector space of dimension d{i). The dimensions of these spaces are explicitly known for low degrees 
(see for example [^]), and we have listed some of them in Table |I]. Finally, notice that, given any 
group G, we have a map 

that associates to every graph F its group theory factor rr{G). This map is an example of a weight 
system for ^(0)'^°'^°. Every gauge group gives a weight system for ^(0)™"^°, but one may in principle 
find weight systems not associated to gauge groups, although so far the only known example is the 
one constructed by Rozansky and Witten in 



which uses instead hyperKahler manifolds. 



i 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


d{£) 


1 


1 


1 


2 


2 


3 


4 


5 


6 


8 



Table 1: Dimensions d{l) of ^(0)^° up to ^ = 10. 
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We can now state very precisely what is the structure of the appearing in 
Feynman diagrams can be grouped into homology classes, we have 



21): since the 



(3.7) 



The factors IriM) appearing in ( |3.7| ) are certain (complicated) integrals of propagators over M. 
It was shown in that these are differentiable invariants of the three-manifold M, and since the 
dependence on the gauge group has been factored out, they only capture topological information 
of M, in contrast to Zk{M), which also depends on the choice of the gauge group. These are the 
universal perturbative invariants defined by Chern-Simons theory. Notice that, at every order i in 
perturbation theory, there are d{i) independent perturbative invariants. Of course, these invariants 
inherit from ^(0)^°°'^ the structure of a finite-dimensional vector space, and it is convenient to pick 
a basis once and for all. Here we will study these invariants up to order 5, and we choose the basis 



presented by Sawon in |39|: 



1: e 



2 : 



3 : 



4 : 



5 : 




(3. 



As in ||3^, we will denote the graphs with k circles joined by lines by Ok- Therefore, the graph 
corresponding to £ = 1 will be denoted by 6, the graph corresponding to £ = 2 will be denoted 62, 
and so on. The second graph for i = 4 will be denoted by u, and the second graph in £ = 5 by 006. 
The group factors associated to these diagrams can be easily computed by using the techniques of 
jT^ (see also @, |3])- Explicit results for all classical gauge groups are presented in the Appendix. 

Remarks: 

1. It is interesting to understand the framing dependence of the universal perturbative invariants 
(see for a discussion of this issue). As shown in the full partition theory Zk{M) changes as 
follows under a change of framing: 

wise 

Z e—Z, (3.9) 
where s G Z labels the choice of framing and 



kd 
k + y 



(3.10) 



7 



is the central charge of the WZW model with group G. Using now that (see Appendix A) 



re{G) = 2yd, (3.11) 

we find that under a change of framing one has 

loiM) ^ leiM) - ^, (3.12) 

while the other universal perturbative invariants remain the same. Since we will work in the 
canonical framing of M, this will produce a canonical value of Ig{M). 

2. Notice that Chern-Simons theory detects the graph homology through the weight system 
associated to Lie algebras. Unfortunately it is known that there is an element of graph homology 



at degree 16 that it is not detected by any weight system associated to simple Lie algebras. However, 
there is a very elegant mathematical definition of the universal perturbative invariant of a three- 
manifold that works directly in the graph homology. This is called the LMO invariant [31| and it 
is a formal linear combination of homology graphs with rational coefficients: 

uj{M)= J2 /r^°(M)r G ^(0)^°"°[Q]. (3.13) 

reyt(0)<;°nn 

It is believed that the universal invariants extracted from Chern-Simons perturbation theory agree 
with the LMO invariant. More precisely, since the LMO invariant uj{M) is taken to be for S^, we 
have: 

as long as the graph F is detected by Lie algebra weight systems. In that sense the LMO invariant 
is more refined than the universal perturbative invariants extracted from Chern-Simons theory. 

3. The Chern-Simons approach to the theory of universal perturbative invariants is very similar 
to the approach to Vassiliev invariants based on the analysis of vevs of Wilson loops in perturbation 
theory |1], |^. The role of graph homology is played there by the homology of chord diagrams (see 
for example [|^, p8|). 

4. The Chern-Simons partition function on Seifert spaces 

In this section we write the partition function of Chern-Simons theory on Seifert homology spheres 
as a sum of integrals over the Cartan subalgebra and a set of residues, by extending results of 
Lawrence and Rozansky ||30| for SU{2). We also show that these integrals can be interpreted in 



terms of matrix integrals associated to a random matrix model. 
4.1 Seifert homology spheres 

Seifert homology spheres can be constructed by performing surgery on a link C in with n + 1 
components, consisting on n parallel and unlinked unknots together with a single unknot whose 
linking number with each of the other n unknots is one. The surgery data are pj/qj for the unlinked 
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unknots, j = 1, ■ ■ ■ , n, and on the final component. Pj is coprime to qj for all j = 1, ■ ■ ■ , n, and the 
Pj's are pairwise coprime. After doing surgery, one obtains the Seifert space M = ■ ■ • , 

This is rational homology sphere whose first homology group Hi{M, Z) has order \H\, where 

n n 

H = pJ2—, and P = JJpj-. (4.1) 

i=i j=i 

Another topological invariant that will enter the computation is the signature of C, which turns out 
to be |3| 

a(/:) = X:sign(|)-sign(|). (4.2) 

For n = 1,2, Seifert homology spheres reduce to lens spaces, and one has that L{p,q) = X{q/p). 
For n = 3, we obtain the Brieskorn homology spheres E(pi,p2;P3) (in this case the manifold is 
independent of qi,q2,(l3)- In particular, E(2,3,5) is the Poincare homology sphere. Finally, the 
Seifert manifold X{^, (^^^^jjyf, ^-y^), with m odd, can be obtained by integer surgery on a (2,m) 
torus knot with framing t. 

4.2 Computation of the partition function 

In order to compute the partition function of M, we first have to compute the invariant of C for 
generic representations P — p, Ai, ■ ■ ■ , A„ of the gauge group G, where P — p is the irreducible 
representation coloring the unknot with surgery data 0, and Aj are irreducible representations 
coloring the unknots with surgery data Pi/qi, i = 1, ■ ■ ■ ,n. This can be easily done by using the 
formula of for connected sums of knots, and one obtains: 

^/3,p+Ai,-,p+A„l'L-j - ■ 14. 

^pP 

Therefore, the partition function of M will be given by 

Z,{M) = e^^- J2 ^^-=^^^+^'^^'_7^+^'^^+^'^ , (4.4) 



where the framing correction is given by the general formula ( |2.9| ). Seifert homology spheres can 
be also obtained by doing surgery on n strands parallel to in x |35|, and then (4^) follows 
from Verlinde's formula [E2 . 



This expression is not suitable for an asymptotic expansion in l/l, since it involves a sum 
over integrable representations that depends itself on /. In order to obtain a useful expression, we 
follow a series of steps generalizing the procedure in [^, First of all, we perform the matrix 

multiplication Ep+A.G^, '5/3p+a,^p+a!'p- ^his gives 

E S*.AA;it^ = exp(=f?s,gu(f (4.5) 
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where the SL(2, Z) transformation in the right hand side is given by 



\ Pi ri J 

and the phase factor is needed in order to keep track of the framing. The partition function is then, 
up to a muhiphcative constant, given by: 

y I 

/ ^ / \n-2 

na>o(^sinf(/3-«)j 

n E E eK)exp|^(-g,/32-2/3(/n, + «;(p))+r,(/n, + Mp))2)|. (4.7) 

If G is simply-laced, the summand is invariant under the simultaneous shift, 

(3 ^ (3 + la, Hi ^ rii- qia, (4.8) 

and also under 

ni^rii+pia. (4.9) 

In these equations, a is any element in the root lattice. This invariance allows us put rij = in the 
above sum by extending the range of /?: [3 = p + la, where p G and a = a^ai, < aj < P. 
It is easy to see that the resulting summand is invariant under the Weyl group W acting on (3, and 
by translations by IPa, where a is any root. We can then sum over Weyl reflections and divide by 
the order of W, denoted by |W|, and use the translation symmetry to extend the sum over /? in 
the above set to a sum over P G (A„//PAr)\A1. Here denotes the set given by the wall of J-'i 
together with its Weyl reflections and translations by IPa inside A^/lPA,. (for SU{N), the wall of 
J-'i is given by the weights with J2iPi — 0- We won't need a precise description of the points of A4 
in the following, since they only enter in the contribution of irreducible flat connections to the path 



integral [^. After performing all these changes, and using the Weyl denominator formula 

a 
2 

a>0 toGW 



Yl 2 sinh - = E e(w)e"'(^), (4.10) 



we can write (|4.7|) as 



1 IZL ^2 V" \ ^ 1 



/3G(A„//PA.)\.M na>o(sinf(/3 



a, 



n-2 



-e-'-ff^' ff r[{-2t) sin -^(/? ■ a). (4.11) 



i=l a>0 
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The last step involves transforming the above sum in a sum over integrals and residues. To do that, 
we generahze slightly and we introduce a holomorphic function of ■ ■ ■ , Pr and xi, ■ ■ ■ ,Xr 
given by: 



Q IP P e ^ ^ 

h{P,x)= ^ — - ^ n~2 Tfr 7i I^2^riaT ^ TT^^ 7i I^2^iwi"' (4-12) 



where j3 = J2l=i Pi^i ^ ® C, x = J2l=i ^i'^i G Ar ® C. This function satisfies: 

h{(3 + IPa, x) = e2'^*^°-^/i(/3, x - IHa), (4.13) 

for any a G A^. Notice also that h{P, x) has poles at the points of A^, the weight lattice. Introduce 
now the integral over C^: 

Q(x) = [ h{P,x)dp (4.14) 

where C" = C x • • • x C is a multiple contour in C, and C is the contour considered in |^: a line 
through the origin from (—1 + i)oo to (1 — i)oo for sign(iJ/P) > (if sign(iJ/P) < 0, we rotate C 
by 7r/2 in the clockwise direction). This contour is chosen to guarantee good convergence properties 
as Pi oo. 

Let us now shift the contour in such a way that it crosses all the poles corresponding to the 
weights in the chamber K^/IPK^. Using ( [4.13| ) it is easy to see that, if Pa ■ x G Z for any root a, 
the resulting integral can be written as 

r r 

^Q{x-lHai)- Q{x-lH{ai + aj)) + --- + {-lY~^e{x-lHj2(^i)- (4-15) 

i=l i<j<i<^ *=1 

The difference between the original integral and the shifted integral ( [4.15|) can be written as 



J2 /(Ax)e-2-*-^. (4.16) 

teAr/HAr ' 



On the other hand, the effect of shifting the contour is to pick the residues corresponding to all the 
weights in the chamber A„//PAj.. Here the residue is understood as lim^-^„- Hjl/^i ~ ni)h{P,x), 
and the residues for the weights that are not in Ai are simply given by {27Ti)~^f{P,x). Putting 
everything together we find, 

J2 f{n,x)= J2 //(A^)e-'^^^'c^/3-(27rO^ J]Res(M/?,x),/3 = n), (4.17) 



whenever Pa ■ a; G Z. We can apply this formula to ( [4.11| ), since what we have there is just a sum 
of expressions of the form /(/5, x) in ( [4.12| ), with x of the form a/P, a G Aj.. In this context, the 
sum over t G Aj./ HA,, is interpreted as a sum over reducible fiat connections on the Seifert sphere, 
and of course t = corresponds to the trivial connection. In the remaining of this paper we will 
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focus on these contributfons, i.e. we will not deal with the residue terms in ( [4.17 ), that should 
give the contribution of irreducible flat connections P^, BDI. In fact, we will only analyze in detail 
the contribution of the trivial connection in order to make contact with the universal perturbative 
invariants. 

In order to present the final result for the contribution of reducible fiat connections to the 
partition function of Chern-Simons theory on Seifert spaces, we have to collect the prefactors, 
including the phases. Define as in |]30| : 



(4.18) 



Therefore, the contribution of reducible fiat connections to the Chern-Simons partition function of 
is given by 



X(Pi 

^91 



:-l)l^+l fVo\A^\ [sign(P)]l^+l^^3,g,(^/p)_^^ 



\W\ (27rz)M VolA, / |P|'-/2 



-e 4 



teAr/HA 



nr=in„>o2sinh|f 



n.>o 2sinh^ 



2p, 
n-2 



(4.19) 



In this equation, (p is given by ( [4.18|) , and in obtaining the phase factor we have made use of the 
Freudenthal-De Vries formula 

1 



p = — dy. 
^ 12 ^ 

We have also introduced the hatted coupling constant 

. Px 



X 



H 



(4.20) 



(4.21) 



where x is the coupling constant given in (|3.3| ). In the evaluation of the above integral we can rotate 
the integration contour C to R*" as long as we are careful with phases in the Gaussian integral, as 
explained for example in If we specialize ( |4.19| ) to G = SU{2), we obtain the result derived in 
||30| . The expression ( [4.19|) is in principle only valid for simply- laced groups, although the results 
for the perturbative series turn out to be valid for any gauge group. 

Notice that, in the sum over A^/HA,., the fs that are related by Weyl transformations corre- 
spond to the same fiat connection. Fortunately, each of the integrals in ( [4.19| ) is invariant under Weyl 
permutations of t, so in order to consider the contribution of a given fiat connection, one can just 
evaluate ( [4.19|) for a particular representative and then multiply by the corresponding degeneracy 
factor {i.e. the number of Weyl-equivalent t configurations giving the same fiat connection). 

If one is just interested in obtaining the contribution of the trivial connection, one can use the 
shorter arguments of and end up with ( |4.19| ) with t = 0. The contribution of the reducible 



connections can also be obtained by generalizing the arguments of [37| to the higher rank situation. 



12 



4.3 Connection to matrix models 

In ( |4.19| ) we have written the contribution of reducible connections to the Chern-Simons partition 
function in terms of an integral over the Cartan subalgebra, since dj3 = YVi=i cind /3j are the 
Dynkin coordinates. In fact, the above expression can be interpreted as the partition function of a 
random matrix model (for a review of random matrices, see ||3^, ^). To see this, let us consider a 
slight generalization of the above results to the U{N) and 0(2r) theories. The partition function 
for these groups can be obtained by writing /3 in terms of the orthonormal basis in the space of 
weights. 

Let us first consider the case of U{N). Denote the orthonormal basis as {ek}k=i, -- ,n, and put 
P = YlkPk^k (where [3k are taken to be independent variables), t = J^k^k^k- It is well-known that 
the positive roots can be written as 

aki = ek-ei, 1 < k < I < N. (4.22) 

Therefore, the integral in (|4.19| ) becomes 

n'' n 2 sinh 
n,<,(2sinh^)" 

We can interpret the /Sk as the eigenvalues of a Hermitian matrix in a Gaussian potential and 
interacting through 

X:$:i«g(2^i-h^^) +(2-n)5^1og(2sinh^^). (4.24) 

i=l k<l ^ A / fc</ ^ ^ 

Notice moreover that for a small separation of the eigenvalues (^4.24 ) becomes, at leading order. 



\og{(3k - Pi f (4.25) 



k<l 



which is the interaction between eigenvalues of the standard Hermitian matrix model. Therefore, 
the integral above can be interpreted as a nonlinear deformation of the usual Gaussian unitary 
ensemble (GUE). In fact, as we will see in detail in the next section, Chern-Simons perturbation 
theory means that we expand around the GUE, and the perturbative corrections are obtained by 
evaluating averages in this ensemble. Note that the non-trivial reducible flat connections, labeled by 
t, are interpreted in the matrix model language as a source term coupling linearly to the eigenvalues. 

Similar considerations apply to the orthogonal group 0(2r). The positive roots can be written 
in terms of an orthonormal basis as follows: 

a^i = ek±ei, 1 < k < I < r, (4.26) 

and the interaction between the eigenvalues reduces again, in the limit of small separation, to 

J2^og{Pl-pfr, (4.27) 

k<l 

which is the eigenvalue interaction of the orthogonal ensemble 0{N) for even N. 
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5. Asymptotic expansion and matrix integrals 

5.1 Asymptotic expansion of the exact result 

In this subsection we will study the asymptotic expansion of the exact result obtained in the previous 
section for the contribution of the trivial connection t = 0. 

The expression (|4.19| ) is very well suited for an asymptotic expansion in powers of x^: we just 
have to expand the integrand in a power series of f3, and integrate the result term by term with the 
Gaussian weight. The integrand has the expansion: 

n„>o 2smh^ a>o 



la>0 \^'^ '^^^^^^ 2 

where /(/?) has the form 

" " " (5.2) 



fi/3) = ll(l + Y,asi(3-ar). 

a>0 ^ s=l ^ 



The coefficients can be obtained in a very straightforward way from (|5.1|) . They are polynomials 
of degree s in n and in the power sums 



Y.P^''- (5.3) 



i=l 



One has, for example, 

1 , 

«i = ^(7ri + 2-n), 

02 = ^—(16 + 5n^ - 18n - lOnvri + 207ri + Svr^ - 2712). (5.4) 
5760 

Let us analyze in more detail the structure of f{f3). Define 

aji(3) = Y,il3-c^y'- (5.5) 

By taking the log of (|5.2| ), one finds: 



(00 
k=i 



00 



4V.(/5) , (5.6) 



where the connected coefficients a^^^ are defined in the usual way: log(l + '^^^k^^) ~ '^k^k^^^- 
An explicit expression for /(/3) can be obtained as follows. Let k = {ki, k2,- ■ ■) be a vector whose 
components are nonnegative integers. Denote i = Yljj^j^ define: 
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Then, 

m = ^+J27-A^^df3^^ (5-8) 

k 

with k\ = Ylj kj\ and the sum is over all vectors k. We see that the perturbative expansion of the 
partition function can be written in terms of the quantities 

^^^^^ - /rf/3A2(/3)e-/^V2 ' ^^-9^ 

where we have denoted 

A\P) = \{{p.af. (5.10) 

Notice that, when we write (3 in terms of the orthogonal basis ( [4.22| ) or ( 4.26| ), ( |5.10| ) is indeed the 
square of the Vandermonde determinant in the variables (3j (for U{N)) or (3j (for 0(2r)). Therefore, 
as we anticipated before, the asymptotic expansion of the integral is an expansion around the 
corresponding Gaussian ensemble, and the perturbative corrections can be evaluated systematically 
as averages in this ensemble. 
We will denote 

' d(3A\(3)e-^"/^, (5.11) 



so that the partition function on Seifert spaces can be written, using ( |4.19|) , as 

log|^ = -^^#x + log(l + f;( J,4''^dG)y)- (5.12) 



=1 



In this equation Zi_ioop is given by 



^i-ioop - |.^| ^2„). y-^J^ j |p|./2 ^ox , (5.13) 

and indeed gives the one-loop contribution around the trivial connection. This follows by comparing 
the exact result with the perturbative expansion 



logf^ = Ef E rr(G)Ir(M))x'. 

^1-loop Vp,^-,„^conn / 



(5.14) 



We also see that, by comparing ( |5.11|) and ( ^.12|) , we can extract the value of the universal pertur- 
bative invariants Ir{M) at each order x^. In order to do that we just have to evaluate TZ^{G) for 

all vectors k with YljJ^j — ^^"^ ^Iso the group factors rr(G) for graphs F with 21 vertices. Of 
course, from a mathematical point of view it is not obvious that the asymptotic expansion of the 
exact partition function has the structure predicted by the perturbation theory analysis. The fact 
that this is the case provides an important consistency check of the procedure. 
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5.2 Evaluating the integrals 

We now address the problem of computing the integrals in ( ^.91) . As we explained in section 4, the 
partition function of Chern-Simons theory on Seifert spaces can be interpreted as a matrix model 
with an interaction between eigenvalues of the form log(sinh(/3j — In the perturbative approach 
we have to expand the sin in power series, and the integrals TZ^{G) are nothing but averages of 
symmetric polynomials in the eigenvalues in a Gaussian matrix model. We will present two methods 
to compute these averages. 

The first method gives the complete answer only up to £ = 5, but it has the advantage of 
providing general expressions for any simply-laced gauge group. The starting point is the following 
identity: 



rf/?e-t^'n4sinh(^^^) sinh 



■ a] 



a>0 

(^)^'V|(det(C))^e%^^^n2-nK^^)' (^-^^^ 



a>0 



where C is the Cartan matrix of the group. This formula is easily proved by using ( |4.1CI| ). Another 



useful fact is that o"i(/?) can be written as (see |]T5l, pp. 519-20) 



Y^{p.af = yP\ (5.16) 



a>0 



One can easily show that, by expanding ( |5.15|) in s, t, and by using ( p.l6| ), it is possible to determine 
the integrals 71^{G) for any gauge group up to ^ = 5, therefore this is enough for the computational 
purposes of the present paper. The answer is given in terms of ?/, ci, and the quantities 

ak = \{{a-pfK (5.17) 

a>0 

For example, one finds: 

7^(o,lA•■.)(^) = 5^^'• (5-18) 
The answers obtained by this method are listed in the Appendix. 

In order to evaluate the integrals ( |5.9| ) for arbitrary o"^, it is important to have a more general 
and systematic method. Here is where the connection to matrix integrals becomes computationally 
useful. It is easy to see that, since the integrals TZj:{G) are normalized, one can evaluate them in 
U{N) and 0(2r) instead of SU{N) and S0{2r). Therefore, one has 

n^^(SUiN)) = dPe-^^f'^/'lliP,- P,Ya^,{P), (5.19) 

i<j 

where 

ani/3) = J2i/3^-/3,)'^. (5.20) 

i<j 
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In (|5.19|) one integrates over N independent variables jSi,--- ,(3^- It is clear that the are 
symmetric polynomials in these variables. One can for example write ( |5.2CI| ) in terms of power 
sum polynomials Pj{P) (defined in ([A.8|) ) as follows, 



1 /9 \ 

s=l ^ ^ 



(/?)• 



(5.21) 



The averages of symmetric polynomials in the Gaussian unitary ensemble can be evaluated in 
principle by using the Selberg integral or the results of [^. A more effective way is the 
following: any symmetric polynomial in the /Sj's can be written as a linear combination of Schur 
polynomials Sx{j3), which are labeled by Young tableaux associated to a partition A (see ( |A.6D ). 
Therefore, if we know how to compute the normalized average of a Schur polynomial. 



df3 e 



n(A-/3i)'^A(/3), 



(5.22) 



we can compute all 7^^. An explicit expression for ( |5.22| ) has been presented in \TA\. The result is 
the following: let |A| be the total number of boxes in the tableau labeled by A, and let Aj denote 
the number of boxes in the i-th row of the Young tableau. Define now the |A| integers as follows 



fi = \i + |A| - i. 



lAI 



(5.23) 



Following we will say that the Young tableau associated to A is even if the number of odd 
/j's is the same as the number of even /j's. Otherwise, we will say that it is odd. If A is odd, the 
normalized average {Sx{j3)) vanishes. Otherwise, it is given by: 



A{A-l) ri/odd/'Tl/'cvcn. -,. , 

-1 2 —J. i— —dim A, 



odd,/' even 



(/ - f ) 



(5.24) 



where A = i/2 (notice that i has to be even in order to have a non vanishing result). Here dim A 
is the dimension of the irreducible representation of SU{N) associated to A, and can be computed 
by using the hook formula. This expression solves the problem of computing the averages ( p.l9|) 
in the general case: we express the product of power sums appearing in (|5.21 ) in terms of Schur 
polynomials by using Frobenius formula ( |A.9| ), and then we compute the averages of these with 
( ^.24| ). As an example of this procedure, let us compute 7^(o,i,o, - )('S'f/(A^)). Using (|5.21| ) and 
Frobenius formula ( |A.9| ), we find: 



a2 = {N- 1)S, 



rrm 



{N + l)Sn + {N -3)S, 



(5.25) 



The averages of the different Schur polynomials can be computed from ( ^.24| ), and we obtain, after 
some simple algebra: 

n^o,i,o,-){SU{N)) = 5N^{N^ - 1), (5.26) 

in agreement with (|5.18 ). 
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Let us now consider the orthogonal ensemble. The averages that we want to compute are given 

by 

^Jd^e-^U^V- n (P'-^lf^m^ (5-27) 



l<i<j<r 



where 



n-1 



= (2r - 22--^)P„(A^) + J2 L)PsiP-)Pn-s{P^). (5.28) 

s=l ^ ^ 

The functions <y^{f3) are now symmetric polynomials in the I3f, so we can write them in terms of 
Schur polynomials S\{(3f). This allows to express the integrals ( |5.27| ) in terms of the integrals 



dy ^\y){y, ■ ■ ■ yrT-'e-^y^^-^y^^/' SM, (5-29) 
which are a special case of a generalization of the Selberg integral studied by Kadell , see also 



33| . Their value is given by 



r! J] T{Xi + a + r-t) J](Ai - A, (5.30) 

i=l i<j 

In our case, a = 1/2. The normalized average of a Schur polynomial is then: 

(SM)) = 2Wdim A n ^'r(V2 + .^\^" (^'^i' 

In this equation, dim A denotes the dimension of the representation of SU{r) associated to A. This 
solves the problem of computing the averages (|5.27| ) in the orthogonal ensemble. As a simple 
example, let us consider again = (0, 1, 0, ■ ■ ■ ). It is easy to see that 

a2 = (2r- 2)5m + (14-2r)5g, (5.32) 

and one finds 

7^(o,l,o,...)('^0(2r)) = 20r(2r - l)(r - 1)^, (5.33) 

in agreement with (|5.18|) . 

5.3 Universal perturbative invariants up to order 5 

Using the above ingredients, it is easy to find the universal perturbative invariants of Seifert spaces 
up to order 5. Although the coefficients in ( |5.2| ) are functions of n and the Newton polynomials 
P^(j9^^), the answer turns out to be more compact when written in terms of elementary symmetric 
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polynomials in the variables ^ — 1 by using (|A.10|) (so for example Ei = —n + Yll=i Pi One 
finds, 

^"'^ 13824 (ff) 

= YI5^^ (;^) (82^4 - 46^3 - 18(1 + E2 + E,)E, - 9Ei - 18E, - QEl - 9j , 

- 663^ (I) ' (^^^^ + ^^^''^^ - ''^ - ^^(^^^ + ''^^ + ^0 ' 
Le = ^^^^[^y {5E, - E,{3E, + 16) + E,{E, + AE, + 12)y (5.34) 

Note that the first universal invariant is given by 

h = , (5.35) 

where A(M) is the Casson invariant of M, in accord with the general result of |36| and with the result 



for the LMO invariant We have also checked that the value for Ig^ listed above agrees with 
the value obtained in § for the LMO invariant using the Aarhus integral. It would be interesting 
to see if these invariants have some nice integrality properties. The perturbative SU (2) invariants 
of integral homology spheres do exhibit some integrality properties (discussed for example in |]30[|), 
but they include extra factors coming from the SU{2) group weights. In general, the invariants 
listed above are rational even for integral homology spheres. For example, Ig.^ G Z/4 for Brieskorn 
integral homology spheres. Also, Ig^ — -^^52 ^ ^/^ those spaces. 

6. Open problems 

Besides the original motivation of understanding universal perturbative invariants from a field theory 
point of view, the results presented here also provide a computationally feasible framework to study 
the Chern-Simons partition function with higher rank gauge groups. There are various avenues to 
explore in the context of Chern-Simons theory and the theory of three- manifold invariants. It would 
be interesting for example to understand the structure of perturbation theory in the background 
of a reducible nontrivial connection, and work out the asymptotic expansion starting from (|4.19| ). 



One could also consider other three-manifolds (not necessarily rational homology spheres) and see 
in particular if the matrix model representation provided here can be generalized to other cases. 
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There are also various interesting physical contexts in which the results of this paper might be 
relevant. Let us end by mentioning a few of them: 



1) The computation of Rozansky-Witten invariants |38| involves the universal perturbative 
invariants of three-manifolds that are extracted from Chern-Simons theory, but the weight system 
is now associated to a hyperKahler manifold (see [^] for a very nice review). Therefore, the universal 



perturbative invariants of Chern-Simons theory are relevant in Rozansky-Witten theory. On the 
other hand, this theory is an essential ingredient in the worldvolume theory of M2 membranes 
in manifolds of G2 holonomy and the computation of the Rozansky-Witten partition function 
should play a role in understanding membrane instanton effects in G2 compactifications of M theory. 

2) Chern-Simons theory on a three-manifold M describes topological A branes wrapping the 
Lagrangian submanifold M in the Calabi-Yau T*M Moreover, the perturbative invariants of 
Chern-Simons theory correspond to topological open string amplitudes on that target. Having a 
systematic procedure to compute Chern-Simons perturbative invariants may prove to be useful in 
further understanding topological open strings in those backgrounds. 

3) Another consequence of our results is that topological A branes in T*M are described by a 
matrix model, when M is a Seifert sphere. It has been recently shown that topological B branes 
on some noncompact Calabi-Yau spaces are described by a Hermitian matrix model characterized 
by a potential iy($) with multiple cuts. It would be interesting to know if there is a relation 
between the matrix models of [1^ and the ones presented here. Notice that according to our results 
the partition function of U{N) Chern-Simons theory on can be written as 

^ e-^^(^^-i) f^dp, -v,.^/2.xTT/. . , 



Jj|^e-S^^?/-n(2sinh^) , (6.1) 



i=l i<j 



This describes open topological A strings on T*S^ with branes wrapping S^, or equivalently 
(after the geometric transition of [^) closed topological strings on the resolved conifold. In the 
limit X — > 0, (|6.1|) gives the standard Gaussian model, as we have argued at length in this paper. 
On the other hand, it is shown in [O that the Gaussian model (which corresponds to Vr(<l') = $^) 



describes type topological B strings in the deformed conifold geometry. This is consistent with 
the fact that, as explained in [^, the deformed conifold geometry gives the mirror of the resolved 
conifold only at small 't Hooft coupling t = Nx, which for fixed A^ means precisely small x. This 
also suggests to consider multicut matrix models with a potential W{^), as in Jl^, but where the 
eigenvalue interaction is not the usual one ni<j(A ~ /^i)^ t)ut n.^j(2 sinh((/5j — /3j)/2))^. In view 
of the above observation, these deformed models might be relevant to understand the mirror of 
the geometric transition studied in |TT], [T^. Indeed, a compact version of this, corresponding to a 
unitary matrix model where the /5j are periodic variables, has been considered in |^ in order to 
describe the stringy realization of the M = 2 Seiberg-Witten geometry. 

4) Although in this paper the focus has been on the perturbative expansion of the partition 
function, the matrix model is also very useful to understand its large A^ expansion. This is an 
interesting problem in itself, and it would be nice to see what is the connection to the approach of 
p!9| . But of course the large A^ expansion of these models is particularly interesting in view of the 
large A^ dualities involving Chern-Simons theory 22|. Although these dualities are not expected 
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to hold for arbitrary Seifert spaces, the results presented here may be useful to understand in detail 
the case of lens spaces (already analyzed in [^) and shed light on the situation for more general 
three-manifolds. 
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A. Appendix 

A.l Group theory factors 

We first present the group theory factors associated to the connected graphs that give a basis of 
^^0^conn order 5. The evaluation of these factors is straightforward by using the graphical 

techniques of Cvitanovic and rather immediate for all of them (except for r^^lG), that gives 



the quartic Casimir in the adjoint and has been computed for all gauge groups in the second 
reference of |]12[)- Our conventions are as in |]12[: the Lie algebra in the defining representation has 



Hermitian generators Ti, i = 1, ■ ■ ■ , d satisfying the commutation relations [Ti,Tj] = iCijk- The 
generators are normalized in such a way that the quadratic Casimir of the adjoint representation 
Ca (which is defined here by C^^ij = X^fc i ^ikiCjki) is twice the dual Coxeter number. This implies 
that Tr(TjTj) = a6ij with a = 1 for SU{N) and Sp{N), and a = 2 for SO{N) (notice that these 
normalizations differ from the ones in 0). 





SU{N) 


SO{N) 


d 


N'-l 


lN{N-l) 


y 


N 


N-2 


Oil 


_LAr2(Ar2_i)(2Ar2_3) 


-^N{N -1){N- 2){8N^ - A5N^ + 5AN + 32) 



Table 2: Dimensions, dual Coxeter numbers and Q2 for SU{N) and SO{N) 



The group factor will be written in terms of the dual Coxeter y, the dimension of the group 
d, and 0^2 (where ak is defined in (|5.17| )). Their values for SU{N), SO{N) are listed in Table 
^. The results for Sp{N) follow from the Sp{N) = SO{—N) relation [|r^, so for Sp{N) one has 
d = N{N + l)/2, y = N + 2 and a2^^^\N) = a2^^^\—N). The group theory factors for the graphs 
in (13.81) are listed in Table |3|. 
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graph 



e 



o — o 
o — o 





group factor 



2dy 



8dy^ 



16dy^ 
I8dy* - 480a2 



32dy^ 
2y{18dy^ - 480a2) 



Table 3: Group theory factors for the Feynman graphs up to ^ = 5. 
A. 2 Matrix integrals 

We now list the results for the matrix integrals (|5.9|), up to order 5. The results for ki = are: 



'^(0,l,0,- 


)(G) 


= bdy\ 


"^(o.o.i.o,-- 


){G) 


= 35 dy^, 


'^(0,2,0,0,-- 


)iG) 


= 25 + 12^ - 288002, 


"^(o.o.oa.o,-- 


)iG) 


= 350(^1/^- 168002, 


"^(0, 1,1,0,0, ■■ 


)iG) 


= 35y{5d{d + 2A)y^ - 1728a2} 


"^(0,0,0,0, 1,0, ■• 


)iG) 


= 4620|/{d|/^- 12a2}. 



(A.l) 



The results for /ci > can be obtained from ( |A.1| ) very easily: insertions of cri(/3) can be reduced 
to insertions of by using (|5.16|) , and these can be computed by taking derivatives with respect 
to a in ( p.l5| ). We have for example: 



7^(2,l,o,...)(G') =5rf(rf + 4)(rf + 6)y^ 
Finally, the integral Zq in ( ^.11 ) is given by 

Zo = (27r)i|>V|(detC)^ H^"''^)- 



(A.2) 



(A.3) 



Q!>0 
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Note that this combines with the rest of the factors in ( |5.13| ) to produce, up to an overall phase 



a>0 



where we have used that VolAj./Vol = det(C). If instead of taking the volume of the root lattice 
in (|A.4| ) we take that of the coroot lattice, the resulting expression for the partition function is 
probably valid for any gauge group. 

A. 3 Symmetric polynomials 

Here we summarize some ingredients of the elementary theory of symmetric functions that are used 
in the paper. A standard reference is [B3 . 



Let Xi, ■ ■ ■ ,xn- denote a set of variables. The elementary symmetric polynomials in these 
variables, Em{x), are defined as: 

Em{x) = Xi,---Xi^. (A.5) 

ii<---<im 

The products of elementary symmetric polynomials provide a basis for the symmetric functions of 
A^ variables with integer coefficients. Another basis is given by the Schur polynomials, Sx{x), which 
are labeled by Young tableaux. A tableau will be denoted here by a partition A = (Ai, A2, ■ ■ ■ , Ap), 
where Aj is the number of boxes of the i-th row of the tableau, and we have Ai > A2 > • ■ ■ > Ap. 
The total number of boxes of a tableau will be denoted by |A| = Aj. The Schur polynomials are 
defined as quotients of determinants, 

deta;^+^-* 

^a(^) = -rrh^- (A-6) 

detXj 

A third set of symmetric functions is given by the Newton polynomials P%{x). These are labeled by 
vectors k = {ki, k2, - ■ ■ , kp), where the kj are nonnegative integers, and they are defined as 

i=i 

where 

N 

P,{x) = Y,<. (A.8) 

i=l 

are power sums. The Newton polynomials are homogeneous of degree ^ = Yljj^j ^"^^ ^ basis 
for the symmetric functions in Xi, ■ ■ ■ ,xn with rational coefficients. They are related to the Schur 
polynomials through the Frobenius formula, 

Pd^) = J2^mS,{x), (A.9) 
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where the sum is over all tableaux such that |A| = £, and Xa(^) is the character of the symmetric 
group in the representation associated to A and evaluated on the conjugacy class associated to 
k (this is the conjugacy class with kj cycles of length j). Finally, one has the following relation 
between elementary symmetric polynomials and Newton polynomials, 

where the sum is over all vectors with ^ . jkj = m. 
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